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Abstract 

We place several new limits on Lorentz violating effects, which can modify particles' dispersion 
relations, by considering the highest energy cosmic rays observed. Since these are hadrons, this 
involves considering the partonic content of such cosmic rays. We get a number of bounds on 
differences in maximum propagation speeds, which are typically bounded at the 10~ 21 level, and 
on momentum dependent dispersion corrections of the form v = 1 ±p 2 /A 2 , which typically bound 
A > 10 21 GeV, well above the Planck scale. For (CPT violating) dispersion correction of the form 
v = 1 -|-p/A, the bounds are up to 15 orders of magnitude beyond the Planck scale. 
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I. INTRODUCTION 



It is generally believed that Lorentz invariance is an exact symmetry of nature. This 
belief is supported by extremely precise experimental tests, and by strong and well motivated 
theoretical prejudice. Indeed, exact Lorentz invariance is used as one of the cornerstones on 
which relativistic quantum field theory is built, leading to the extremely successful Standard 
Model of Particle Physics. Together with the equivalence principle, local Lorentz invariance 
is one of the assumptions underpinning General Relativity. Alternatively, we could say that 
General Relativity is the gauge theory of general coordinate invariance, and has exact local 
Lorentz invariance as a consequence. In either case, it is intimately related to both our best 
theories of particle physics and of gravity. 

While most physicists believe that Lorentz symmetry is an exact symmetry of nature, one 
of our jobs as physicists is to put all of our assumptions to the most vigorous possible tests. 
Therefore, we should attempt to place the tightest possible limits on the violation of Lorentz 
symmetry. A general framework for considering Lorentz violation has been developed by 
Kostelecky and collaborators 0. The best laboratory experimental limits come from clock 
comparison experiments, for instance the recent experiments of Cane et. al. (H, see also 
These experiments are highly constraining, particularly for low energy effects and those 
involving CPT symmetry violation. (In quantum field theory, CPT symmetry is a conse- 
quence of Lorentz symmetry, together with causality, localitvanalyticity, and unitarity, and 
its violation implies that Lorentz symmetry is also violated |4|.) 

However, it is consistent and reasonable to believe that CPT could be an exact symmetry 
without Lorentz invariance. It is also desirable to obtain limits involving rather high energies. 
In fact, it turns out high energies introduce another way of gaining high precision, in a 
manner suggested by Coleman and Glashow p, |6| . 

The argument goes as follows. If Lorentz symmetry is not exact, then different particles 
can have different limiting velocities, and potentially more complicated corrections to their 
dispersion relations. This opens up the possibility of Cherenkov radiation; a particle with 
a larger limiting velocity is energetically allowed to radiate particles with lower limiting 
velocities, even though such radiation by a particle at rest is not kinematically allowed. For 
instance, suppose the limiting velocity of an electron is c e > c 7 , where c 7 is the limiting 
velocity of a photon. Then the energy of an electron of momentum p ^> mc e is 

2 3 

E e , p = ^cl + mlci ~pc e + ^, (1.1) 

Zp 

while the energy of a photon of momentum k is kc-y. The process e — > is kinematically 
allowed if the final state energies sum to the initial state energy and the sum of the mag- 
nitudes of the final state momenta exceeds the magnitude of initial state momentum (so a 
nonvanishing opening angle is allowed). Denoting the final photon momentum as k, defining 
e = (c e — c 7 )/c 7 , and taking m 2 e c 2 e /p 2 ~ e < 1, this condition becomes 

Cherenkov condition: ke > — ^— for some k, or e > ^- , (1.2) 

2p{p — k) 2p 2 

which is the same as the condition that the electron is traveling faster than the speed of 
light c^. This is "traditional" Cherenkov radiation. 
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The highest energy particles observed or inferred are not those accelerated in laboratory, 
but those which occur in cosmic rays or other astrophysical settings. The Cherenkov process 
is so efficient that an electron which exceeds the threshold for Cherenkov radiation would 
lose energy down to the threshold on a time scale which is very short compared to any astro- 
physical time scale. Therefore, electrons of higher energy would never exist astrophysically. 
One can then place a limit on e of 

mlc 2 

e<^, (1.3) 

with p the momentum of the most energetic electron observed or inferred in an astrophysical 
setting. (Since the constraint on e will be so tight, there is no need to specify which c appears 
on the right hand side.) 

Similarly, if the photon's limiting velocity is larger, then the process 7 — > e + e~ becomes 
kinematically allowed for photons of momentum k satisfying 



Pair production condition: /cc 7 > kc e H f — , (1.4) 



27BgC 3 

T 



where we already used that the process is kinematically most favored when the electron and 
positron share the momentum evenly. This leads to a bound, 

with k the momentum of the most energetic photon observed. The process involved here is 
not "traditional" Cherenkov radiation, but a sort of "generalized" Cherenkov process. 

Observations of cosmic ray photons, and of photons inferred to arise from synchrotron 
emission off electrons in the Crab nebula, give the constraints 0| 

-2 x 1(T 16 < e < 5 x 1(T 20 , (1.6) 

a two-sided bound on Lorentz violation (though not as strong as those from clock experi- 
ments). The authors of j?} also show that such Cherenkov limits place bounds on dispersion 
corrections involving additional powers of the particle momentum, of form E = pc + p 2 /A, 
which limit A to be at or above the Planck scale (a constraint much stronger than can be 
obtained from existing laboratory experiments). 

One would expect that the strongest limits of this sort would arise from the cosmic rays 
with the largest energy. The most energetic cosmic rays observed carry energy in excess of 
10 11 GeV and are believed to arise from hadronic primaries 0. The highest energy cosmic 
ray observed had an energy of about 3 x 10 11 GeV [9], and was probably a hadron (lp| . 
The fact that such high energy particles can arrive from astronomical distances gives the 
constraints j^i HI)]. 

C p — C 7 23 



Cp 



< 10~ 2 % (1.7) 

< 1(T 15 , (1.8) 

Cp 

where c grav is the propagation speed of gravitation. (The constraint involving gravitation is 
weaker because gravitational Cherenkov radiation is inefficient; the velocity difference must 
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be 10~ 15 before gravi-Cherenkov radiation becomes appriciable on an astrophysical time 
scale.) Unfortunately, these bounds are both one-sided; c p could be substantially slower 
than c 7 or c grav without violating either bound. 

The proton is a composite particle. At the energy relevant to the current problem, 10 11 
GeV, the proton (as viewed in our rest frame) can be considered as a loose bag of several 
partons, mostly quarks and gluons, each carrying a small fraction of the proton's total energy. 
The purpose of this paper is to show that, by taking into account the partonic structure of 
hadrons, we can use the arrival of high energy hadrons from astronomical distances to place 
a number of severe limits on Lorentz violating physics. For the reader's convenience, we will 
present the main results here. 

In all of our results, we will assume that the highest energy extensive air showers arise 
when single hadrons strike the upper atmosphere. (It is possible that heavier nuclei are 
actually responsible; this would weaken some of our bounds.) For some results, we will 
make an additional assumption; in a "bottom-up" generation mechanism for the high energy 
cosmic rays. That is, we will sometimes assume that the highest energy cosmic rays got their 
energy by being accelerated by ordinary electromagnetic fields over astronomical distances, 
rather than, for instance, arising as the decay products of some long lived and ultra-heavy 
relic particles. In this case, the hadron must have been electrically charged during the 
acceleration process. Therefore, some charged hadron must be stable, at least on the time 
scale of hours, at the energy of the highest energy cosmic rays observed. In what follows, 
we will label results which depend on this assumption with an asterisk *. 

Differences in propagation speeds arise from dimension 4 Lorentz violating operators. We 
will parameterize these by the propagation speed difference. Choosing to renormalize such 
operators at the scale 10 11 GeV (see below), we find the following limits: 

-1.6xl(T 23 <^i^< 1.4 xlO- 23 , (1.9) 
* - 1.8 x 1(T 21 < £2 ^ < 1.2 xlO -21 , * (1.10) 

* - 2 x 1(T 21 < ( c «+ c 3)/2-c 7 < 5 x 1Q -24 ^ ( 1U ) 

* - 1 x 1(T 20 < (c * + ^ /2 ~ Ce < 5 x 1(T 24 . (1.12) 

In the first line, c\ and Ci/ 2 are the propagation speeds of fundamental gauge bosons and 
fundamental fermions, if propagation speeds are common to particles of a given spin but 
differ for particles of different spins. In the second line, c q is the propagation speed for quarks 
assuming all quarks share a maximum velocity, and c g is for gluons. In the remaining lines, 
c 7 and c e are the propagation speeds for photons and electrons (the physical particles) 
respectively; the limits on electrons also apply for neutrinos. 

Note that all of these limits are two sided, and many do not require an assumption 
about the nature of cosmic ray acceleration. These constitute severe constraints on Lorentz 
violation. In deriving these limits we have assumed that the 3 x 10 11 GeV cosmic ray was a 
hadron and that its energy was correctly measured; the strength of the bound involves the 
square of the cosmic ray's energy, so if the highest energy cosmic ray hadron energy were 
only 1 x 10 11 GeV, the results are an order of magnitude weaker. 

Dimension 5 operators induce spin dependent 0(p 2 /A) corrections to particle dispersion 
relations 0| , with opposite sign for particle and anti-particle (and for opposite polarization 
states of gauge fields). There is an independent such parameter for each Standard Model 
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field. Note that these operators violate CPT, and so can be avoided by insisting on CPT 
symmetry. 

Unless the coefficient k of these corrections, E = p ± kp 2 , are at least 20 times larger 
for QCD degrees of freedom than for photons and leptons, we can obtain the following 
constraints on k for photons and leptons: 
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These limits are about 15 orders of magnitude stronger than the Planck scale, and do not 
depend on assumptions about the origins of cosmic rays. This means that any theory 
which induces dimension 5, CPT violating operators must generate them at a highly trans- 
Planckian scale, or must introduce them with a very small coefficient ^ 10 -15 . However, for 
certain peculiar corners of the parameter space of dimension 5 operators, we cannot place 
limits. For instance, if the right-handed up quark receives the largest dispersion correction 
and it is in the range E = p — k UR p 2 , with k > 1 /(10 32 GeV); and if dispersion corrections of 
all non-QCD particles are at least 30 times smaller, then the A ++ may be the most stable 
propagating particle at the energies of the highest energy cosmic rays, and no limits can be 
placed. 

We are able to place super-Planck scale limits (assuming the Lorentz violating physics 
introduces the operators with 0(1) coefficients) even on dimension 6 operators. Such oper- 
ators induce CPT-even dispersion corrections of form E = p + Kp 3 . The limits we obtain 
on the coefficient K are, 
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or 25 1 lf g | *, 

or 500(^1 * . (1.14) 

Here K x and K 1 / 2 are the dispersion corrections of spin-1 and spin-1/2 particles if these 
are taken to be common to all particles of the same spin; K q is the dispersion correction 
of quarks if it is taken to be flavor independent, and K\ cpt is the largest of the dispersion 
corrections of any light lepton (e, /i, z/ e , z/ M , or u T ). In the last two lines, the bound is 
whichever limit is weaker. 
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We believe that these constraints will make it next to impossible to develop a theory in 
which Lorentz symmetry is entirely absent, since they imply a scale of violation at least 
2 orders of magnitude higher than the Planck scale. They are rather less constraining 
if Lorentz symmetry is spontaneously broken in some hidden sector. Note that the scale 
appearing in the bound arises as p 2 /m p , with p the momentum of the highest energy cosmic 
ray; the number quoted should be scaled as the second power of that energy, if 3 x 10 11 GeV 
does not hold up as the energy of the highest energy cosmic ray. Some of these bounds, 
such as the lower bounds on K\ and K 1 , apply even if the cosmic rays are nuclei rather than 
individual hadrons. 

The remainder of the paper is organized as follows. In Section |HJ we review and discuss 
what Lorentz violating operators look like and how they generate corrections to dispersion 
relations. In Section IIII1 we show that the corrections to a hadron's dispersion relation 
involve the expectation value of the operator in a state given by a highly boosted hadron. 
This in turn is related to a twist two operator which has a simple expression in terms of 
the parton distribution functions (PDFs) of the hadron, evaluated at a scale set by the 
parton's energy in our frame. In Section IIV[ we present our results for integral moments of 
parton distribution functions renormalized to the scale 10 11 GeV. Section fVl shows how the 
resulting PDFs can lead to specific decay processes which generate the constraints we have 
just presented. We end with a short conclusion. Finally, details on PDF evolution in the 
Standard Model are postponed to an appendix. 



II. LORENTZ VIOLATING OPERATORS 



It is fair to assume that any Lorentz violating effects are small, in which case we can 
treat Lorentz violation as a small perturbation on Lorentz invariant ordinary physics. In 
the language of field theory, Lorentz violation can be described by the introduction of Lorentz 
violating operators into the (Standard Model) Lagrangian. 

Fairly general discussions of Lorentz violating operators have been presented by Colladay 
and Kostelecky in a general context, and by Coleman and Glashow in the context of the 
Standard Model (if. The most convenient way to introduce such operators is to introduce 
one or more preferred frame 4-vectors. One then constructs Lagrangian terms as usual, 
except that Lorentz indices can be contracted against preferred frame 4-vectors as well as 
against each other. 

Consider for instance electrodynamics with a fermionic field if) and a gauge boson 
with a single timelike 4-vector w M , u^u^ = —1, defining a preferred frame 1 . Besides the usual 
Lagrangian terms, 

rmpip + tp^D^tp + -^F^F^ , (2.1) 
we can write the new dimension 4 terms, 

c" 2 - 1 

(c" 1 - Vju^u^D^ + > ( 2 - 2 ) 



We use [ — h + + ] metric convention. 
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which modify the free velocities of propagation, in the frame where is a unit time vector, 
as indicated. We can eliminate one or the other term by performing a rescaling of the 
coordinate system, stretching the "time" coordinate (the coordinate which is time in the 
frame where vP 1 points in the time direction), so that one particle velocity again becomes 
1. However, the ratio of velocities for the different fields does not change under such a 
coordinate redefinition, so if this ratio is not 1, then Lorentz violation cannot be eliminated 
by any coordinate transformation. 2 

It is also possible to write Lorentz violating operators of both lower and higher dimension. 
For instance, a dimension 3 term, 

tVM^, (2.3) 

splits the energies of the particle and antiparticle, a CPT violating effect. At dimension 6, a 
term of form —Ku tJi u a upu 1 il) r j^d a d l3 d 1 ip modifies the dispersion relation, so that E — KE 3 , 
rather than E, appears (in the frame where w M is purely temporal). At high energy and 
treating the perturbation to be small, the dispersion relation in the rest frame is then 
approximately E ~ p + Kp 3 , which perturbs the propagation velocity by v — 1 + 3Kp 2 
(setting c = 1). 

The coefficients of Lorentz violating operators will run with scale, and will generically 
mix. Rather than investigate this in detail, we will constrain their values at the scale set by 
the physics which is providing the constraint, which in our case will be ~ 10 11 GeV. 

The absolutely most general Lorentz violating Lagrangian has enough terms to make the 
analysis very complicated. We will make a number of simplifying assumptions: 

1. There is a unique timelike 4- vector appearing in Lorentz violating Lagrangian 
terms, and its preferred frame is not too different from the frame of the microwave 
background (the boost factor is at most 0(1)); 

2. For the most part, we will only consider terms even in m m . Terms odd in violate 
CPT, giving opposite corrections to a particle and its CP conjugate. 

3. Lorentz violating effects will be family diagonal. We will also present certain con- 
straints assuming they are family blind, or even that they are the same for all gauge 
bosons and for all fermions. 

The last assumption is for convenience, and because we don't know of any strong motivation 
for theories where it is violated. The other assumptions do not seriously restrict our con- 
clusions, because the cases where they do not hold are typically more strongly constrained. 
For instance, dropping CPT symmetry allows dimension 3 operators, on which there are 
very tight laboratory constraints 0]. Nevertheless, we will consider what happens when it is 
violated-though our ignorance of spin dependent parton distribution functions will make it 
hard to make as complete an analysis as we make for the CPT symmetric case. Regarding 
the first assumption, we make two comments. Constraints on spacelike 4-vectors are about 



2 It is worth remarking that other considerations, such as causality and locality, may limit what terms 
can appear in such a Lagrangian [lflj . For instance, Kostelecky and Lehnert argue that the term we have 
written, with v < 1, cannot be canonically quantized in all frames, while the v > 1 case may have causality 
violation issues 
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as strong as those on timelike 4-vectors because we have observed ultra-high energy cosmic 
rays arriving over a wide swath of the sky. For the case where is purely temporal (or 
purely spatial) only in a frame of reference at a large boost with respect to our frame, some 
of the cosmic rays we observe have much higher energies in that frame, and would therefore 
Cherenkov radiate even more strongly. Therefore such are more tightly constrained than 
the u^ 1 which is at rest in the microwave background frame. 

To summarize, Lorentz violating effects behave, at dimension 4, like species dependent 
speeds of light, and at dimension 6, like shifts in the propagation speed of form 8v ~ p 2 /A 2 . 
Because it is convenient and intuitively clear, we will generally parameterize the Lorentz 
violating operators in terms of their dispersion relation corrections in the remainder of the 
paper. 



III. HADRON DISPERSION AND THE TWIST EXPANSION 

How does a Lorentz violating perturbation to the Lagrangian modify a particle's 
dispersion relation? 

First, we must derive from the perturbation to the 4- momentum operator P^. The 
perturbation to the particle's 4-momentum is then given by 

<5p" = MP|^llM, (3-1) 

where tp is the particle type and p is the particle's momentum. The velocity is v = dp°/d\p\, 
so it is modified, in the ultra-relativistic and small Lorentz violation limit, by 5v ~ d(5p° — 
5|p|)/d|p|. Since we can shift by a Lorentz conserving piece, we can generally make 
S\p\ = 0, in which case the correction is just dSp°/d\p\. Evaluating this for a free particle 
is straightforward; we find the part of which looks like a kinetic term for this particle, 
and replace derivatives with the momentum p^. The free particle relation between the 
Lagrangian terms and the modification of particle dispersion relations then holds. 

But what if the particle in question is a hadron? Naively, we should view the hadron 
as a sum of partons with different momentum fractions x of the hadron's momentum, and 
expect 8p^ to be the sum of the contribution expected from each parton. That is, suppose 
that the dispersion relation for parton type / is 

Sp°[fM = </,p|^° l/,P> = (c f - l)|p| +^/|p| 3 . (3.2) 

Write the parton distribution function for a hadron h, such that the probability to find a 
parton / with momentum fraction between x and x+dx when using a probe of virtuality q 2 
is given by fh(x, q 2 )dx. We should then expect the correction to the hadron's energy to be, 

5p°[h,p] = (h,p\P^\h,p) = J2 [ dx\(c f -l)x\p\+K f x 3 \p\ 3 ]f h (x,\p\ 2 ). (3.3) 

/ J 

This basically turns out to be the case. 

For our case, Py is simplest in the rest frame of the 4- vector u^, and in this frame, |p| 3> 
Aqcdj showing that the expectation value of the Lorentz violating operator is determined 
by hard physics. Alternatively, we could say that u^p^ provides a hard scale which makes 
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a partonic description applicable. When evaluating an operator such as in the wave 
function of a boosted hadron, as in Eq. I)3.ip . one should perform an operator product 
expansion on in terms of low twist operators. The dominant contribution will arise from 
the twist-two part of the operator, with higher twists suppressed by powers of Aqcd /u^p^ ~ 
10 |; . 

Consider a twist two, dimension D operator Od- Denote its expectation value in an 
unphysical state with parton type / carrying momentum p to be 

(f,p\0 D \f,p) = (On)f. (3.4) 

Then the definition of the parton distribution functions (PDFs) is, that the expectation 
value of in a hadron is given at leading order by 

(h,p\0 D \h,p) = ^2(0 D ) f [ dxx D ' 3 f h (x) + 0(a s /2n), (3.5) 

where fh(x) is the PDF for parton / in the hadron h. This is valid up to higher order effects, 
which are minimized (order a s /27c without logarithmic enhancement) if we take the operator 
O to be renormalized at the relevant hard scale, which for us is about 10 11 GeV. 3 In other 
words, the expectation values of twist-two operators are given by the Mellin moments of the 
PDFs-or to be more accurate, at leading order the PDFs are defined as the inverse Mellin 
transforms of the twist-two moments of the hadron. 

Therefore, the correct statement is that the value of an operator in a hadron is the 
sum over partonic species of the integral over x of the PDF times the expectation of the 
operator for the given parton. In other words, our naive expectation in Eq. ()3.3|) is correct, 
up to 0(a s /27r) corrections. Since ot s ~ 0.03 at the scale in question, we can safely ignore 
such corrections. (The corrections could be found by making a loop-level treatment of the 
operator product expansion of the Lorentz violating operator in the twist expansion, and a 
loop level evaluation of the twist-two operator in terms of the PDFs of the hadron.) 

Therefore, the effect of dimension 4 Lorentz violating operators is to give a hadron a 
limiting velocity of propagation of 

Ch = c f I dxxfh(x), (3.6) 

with the sum running over all partonic species. (Note that Jq dx xfh(x) = 1 to ensure 
momentum conservation.) The limiting velocity of a hadron is the momentum weighted 
average of the naive limiting velocities of its partonic constituents. Similarly, if parton type 



3 The Lorentz violating operators are presumably induced at a scale of order the Planck scale and must 
then be RG evolved from that scale to the scale 10 11 GeV. Such evolution has been studied for the case 
of QED in The evolution generically leads to 0(a log(A//x)) mixing between operators, making it 
unlikely that one operator should be smaller than another of the same dimension by more than about a 
factor of 10. One also expects high dimension operators of dimension D to induce lower dimension ones of 
dimension d with coefficient ~ aA D ~ d . This is suggests that dimension 4 operators will have coefficients 
0(a) 15], which as we will see is grossly in conflict with observation. 
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/ naively gives a dispersion relation of E = p + Kjp 3 due to dimension 6 Lorentz violating 
operators, then the energy of a hadron will be given by 

E = p + p>Y^K f [ dxx 3 f h (x). (3.7) 
/ Jo 

Since most of the momentum of the hadron is carried at relatively small x, this means that 
the correction to the hadron dispersion, due to dimension 6 operators, will be much smaller 
than it would be were the hadron not composite. 



IV. STANDARD MODEL PARTONIC CONTENT OF HADRONS 



The evolution kernel for parton distribution functions, within QCD, is known to next-to- 
leading order (NLO) |1(tI| , and the evolution equations for lower order integral moments are 
known to NNLO [13, lla]. Global fits of the PDFs of the proton have been performed using 
data at many energy scales and values of x [l^, HtJ 21 1; those of the neutron are given, to 
suitable accuracy, by replacing u d. None of the available PDF libraries go to 10 11 GeV, 
and none include the evolution of electroweak degrees of freedom (the PDFs of photons and 
weak bosons, leptons, and the Higgs boson within the proton). 

Since some of our most interesting results will come from analyzing the relatively small 
contamination of electroweak degrees of freedom in hadrons, it behooves us to extend the 
analysis in this way. Luckily, we are only interested in integral moments of the PDFs, 
which actually have much simpler evolution equations with scale. Namely, the integral 
moments of form J Q x n f(x)dx obey a closed set of ordinary differential (renormalization 
group) equations, which can be evolved along with the renormalization group (RG) equations 
for the couplings, once initial values for the integral moments are determined from the global 
fits. The appendix explains why the integral moments satisfy ordinary differential equations 
(a point which is well known at least to parts of the community), and extends the usual 
treatment to include the LO effects of all Standard Model interactions, including electroweak 
interactions with scalars and Yukawa interactions. 

Since we are not generally interested in tracking how much of the structure is in particles 
and how much in ant-particles, we will only keep track of the sum of the PDFs, q(x) + q(x) 
(the singlet distribution), which we will nevertheless label q(x). We will also not track spin 
dependent PDFs. The spin dependence and particle-antiparticle difference are both needed 
to place detailed limits on dimension 5 CPT violating operators, but the spin dependent 
PDFs are not known with enough accuracy to justify such a treatment anyway. 

We have evaluated the relevant integral moments for the quark and gluon PDFs in a 
proton, neutron, andpion at the relatively low scale of 3 GeV, by using a fit from the 
CTEQ collaboration 3| for the p, n and an old MRST collaboration fit for the pion. 
These were evolved to the scale 10 11 GeV. In this evolution, gauge field beta functions and 
QCD contributions to the DGLAP equations (the equations describing the scale dependence 
of the PDFs) were treated at the two loop level, while the running of the Yukawa coupling 
and all other DGLAP evolution was treated at one loop level. The Higgs self-coupling was 
not treated at all, because it first leads to DGLAP evolution at two loops (and then only 
in the uninteresting Higgs scalar sector) and to Yukawa and gauge boson beta functions 
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TABLE I: Dimension 4 content, at 10 11 GeV scale, of p, n, 7r, 7, e, D°, and glueball (gg). 

at two and three loops respectively. We also treated all thresholds in a rather naive way; 
for mass thresholds, a particle was excluded from the theory until q 2 = m 2 , and then 
included as massless. For the electroweak threshold, taken to occur at mz = 91.19 GeV, 
we abruptly switched to treating QED+QCD to treating the Standard Model, neglecting 
all particle masses. At this threshold, the PDF for the u quark was split evenly between 
the right and left handed components (appropriate for a spin averaged sample of protons), 
and then the left u type and d type quarks were mixed with each other, while the photonic 
PDF content was mixed into the B and W boson types according to cos 2 Qw and sin 2 Q\y 
respectively. (Mixing left quarks would occur anyway due to soft W radiation, a process 
which is logarithmically divergently efficient.) 

This treatment of PDF evolution, as well as the starting PDFs, naturally introduce errors. 
This is particularly true for the rather poorly known pion PDFs. Note however that the 
PDF evolution of the integral moments in question tends to converge towards a fixed point, 
particularly for the strongly interacting degrees of freedom (because of the large size of 
a s and the consequent rapid evolution of the PDFs). This is evidenced by the startling 
similarity of the proton and pion PDFs at the ultraviolet scale. As a consequence, the errors 
are reduced. In addition, in what follows we will be most interested in the difference between 
the PDFs of different particles, and the errors in these differences-particularly differences 
arising because of electroweak interactions-are much smaller. For these reasons, we will 
not make a detailed error analysis of our PDF evolution. (In any case, our final errors will 
be dominated by the error in the energy of the highest energy cosmic ray observed to date, 
which may introduce as much as a factor of 2 error in the numbers in our bounds.) That said, 
we should particularly point out that the pion PDFs are not as well known, in particular 
the size of the quark sea is uncertain. Using the extreme large and small sea values given 
in ji^l, we find that the gluon fraction in Table |l] varies by ±.004. Fortunately, we do not 
find that processes involving pions are particularly important in bounding Lorentz violating 
operators. 

We also evolved the PDFs for the electron, the photon, and crude models of the D meson 
and the glueball. For the electron and the photon, we took the partonic content to be a 
pure electron or photon with x — 1 at the scale 500 KeV, and evolved under electromagnetic 
interactions only to the scale 500 MeV, and with full interactions from there. This procedure 
will underestimate the QCD content of the photon somewhat; but we find that neither the 
photon nor the electron takes on much admixture of other degrees of freedom, and it is not 
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TABLE II: Dimension 6 content, at 10 11 GeV scale, of p, n, ir, 7, e, D°, and glueball (gg); in 
percent (%) of value for a particle with the full x = 1 fraction of the momentum. 

important in practice to use the evolved rather than pure structure. 

Treating the D and the glueball requires more modeling. We treat the D meson, at 
the scale m c ~ 1.5 GeV, as a charm quark with x = 0.8, with the remaining 20% of the 
momentum made up equally of light quark and gluon. The light quark and gluon content 
are taken to have {x 2 } = 0.01 each, so they contribute quite little to the high dimension 
content of the particle. The structure at higher scales is then obtained by DGLAP running, 
starting at the scale 1.5 GeV. We considered similar models of the J/ip and B mesons, but 
these turn out to be too heavy to be important in setting the Lorentz violating bounds. The 
glueball was taken to be pure glue at the scale 0.5 GeV, with (x 2 ) = 0.2, and was evolved 
from there. We think this model of the D° is reasonable, but the glueball model is more 
speculative. 

The integral moments we consider can in principle be determined by lattice gauge theory 
techniques. This would be particularly useful for the glueball, as otherwise we have little 
guidance in how to treat it. 

We find that all QCD particles are about 45% glue and 50% quarks, with the remainder 
made up mostly by gauge bosons. The glueball is richer in glue and the D richer in quarks, 
but not strikingly so. The photon and electron are about 80% the original constituent. As 
for the x 2 weighted structure, in every case the sum over constituents gives less than 1 (as 
it must; splitting conserves ^2x but always reduces (a; 2 ).) However, for the electron and 
photon, it is a modest reduction, while for hadrons the reduction is enormous; weighting with 
x 2 , barely anything is left (or alternately, most of the partonic content has x < 0.1). This is 
particularly so for the glueball; the root-mean-squared value of x for a glueball constituent 
is about 0.03. This is because gluon radiation, which brings down a particle's x, is efficient, 
especially for gluons themselves. 
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V. DERIVATION OF MAIN RESULTS 



A. Dimension 4 operators 



Consider first, constraints on dimension 4 Lorentz violating operators. These look like 
differences, between particles, in the maximum speed of propagation. Since our constraints 
will always require the Lorentz violating operators to be subdominant to the Lorentz con- 
serving ones, and since the particles involved will always be highly relativistic, we will always 
work to linear order in velocity differences, and will not distinguish different particle veloc- 
ities from each other, and will set c = 1, after we have made the leading order expansion. 
Similarly, masses will be taken to correct dispersion relations only at leading order in m 2 /p 2 . 
So a particle with propagation speed (1 + e) and mass m will be taken to have energy 

m 2 



E = y/(l + 7) V + (1 + e) 4 m 2 = p + tp + — . (5.1) 

2p 

Consider first a difference in velocity between fermions and gauge bosons; all fermions 
have a common propagation velocity at the scale 10 11 GeV, and all gauge bosons have a 
different common velocity. This could happen if the source of Lorentz violation, whatever it 
were, distinguished spin but didn't "care" about gauge charges or other internal symmetries. 

In this case, it is quite easy to find a two sided bound on the propagation speed difference. 
Call that difference, e = c\ — C\i% (with c\ the gauge boson speed and Cu 2 the fermion speed). 
Consider first the case e < 0. In this case, the process p — > pj is kinematically allowed at 
suitably high energy. The proton content is almost exactly half each, gauge boson and 
fermion, whereas the photon is about 80% gauge boson and only 20% fermion (see Table 
P). The energy difference between a proton of momentum p, and a final state with collinear 
proton of momentum p(l — x) and photon of momentum px, is 

2 2 
771 771 

AE = (1 + 0.5e)p + — E - (1 + 0.5e)(l - x)p , p ? - x(l + 0.8e)p 

2p 2p(l — x) 



P 



-0.3ex 



x m 2 



x 2p 2 



(5.2) 



When this difference becomes positive, the initial state energy is larger than the final state 
energy, meaning that the decay can occur with a finite opening angle (which will make the 
final state energy a bit larger and use up the extra energy). The optimal value of x for this 
to happen is x ~ 1, and the process becomes kinematically allowed at a threshold, 

mt 

e = V (5.3) 

When this bound is exceeded by a factor of 2, the proton loses of order half its energy on 
a length scale of less than a meter; so the bound on cosmic rays can be taken essentially to 
coincide with the kinematic bound for the process. 

Now suppose e > 0, so fundamental gauge bosons propagate faster than fundamental 
fermions. In this case, contrary to simpleminded intuition, the breakup of a proton, p — > 
pe + e~ , becomes allowed, since the electron is about 85% fermion and only 15% gauge boson 
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(almost exclusively electroweak gauge boson). The kinematics are almost the same as for 
the previous case, since the electron mass is negligibly small. Combining the two, the bound 
on e is, 

-^<(ci-c 1/2 )<^. (5.4) 
0.6p 2 V 1 ' 0.7p 2 V ' 

Substituting in m p = 0.94 GeV and p = 3 x 10 11 GeV gives the bounds quoted in the 
introduction. Note that neither bound relies on any details of the cosmic ray production 
mechanism; all we need to know is that the highest energy cosmic rays are hadrons, and 
that all hadrons are essentially even mixes of gauge boson and fermion. 

We have not worried, in the above limit, about scalar (Higgs boson) content, because 
no low energy particle has much Higgs content. If for some reason c\ and Ci/ 2 were ap- 
proximately equal and cq, the fundamental scalar propagation speed, were more different, 
we could also achieve limits, about 200 times weaker, by the same processes, because the 
photon contains more, and the electron less, Higgs boson than any hadron made of light 
constituents. 

These limits can be evaded if we assume that propagation speed differences are not uni- 
versal but show more complicated species dependence. For instance, suppose that strongly 
interacting particles propagated slightly more slowly than other degrees of freedom. We can 
still place a bound, albeit a weaker one, on the propagation speed difference between quarks 
and gluons, e' = (c q — c g ), as follows. 

Suppose first that gluons propagate more slowly. Then a gluon-rich hadron would be 
energetically more favorable. A hadron of mass m\ and gluon fraction g\ would split off a 
daughter hadron of mass m 2 and gluon fraction g 2 provided the energy difference, 

AE = (1 - 9l e')p + ^ - (1 - gxe')(l -x)p- n f 1 . - (1 - g 2 e')xp 



2p ' 2p(l — x) 2px 

2 2 

= x{92 - 9l)ep -^p-T^-2p^ (5 - 5) 

is positive for any x. The most gluon-rich hadron that we have been able to identify is 
the glueball, a ++ scalar. There is no conservation rule to forbid the process p — > p + gg 
(gg the glueball). Taking the glueball mass to be 1.6 GeV [23], we find that the process 
goes forward, for a proton with energy equal to the highest energy cosmic ray observed, if 
e' > 1.2 x HT 21 . 

The presence of the process p — > p + gg does not necessarily mean that cosmic rays could 
not reach the Earth; they would simply have to be composed of glueballs. Even though 
glueballs are unstable at rest, in the Lorentz violating scenario under discussion, they would 
be stabilized by the Lorentz violating terms when their energy reached that of the highest 
energy cosmic rays. Therefore they could propagate to the Earth over cosmological distances. 
This is provided they were not destroyed by the microwave background; but in fact this seems 
likely, because a glueball should have a small cross-section with photons. Therefore, this 
scenario would probably even evade the GZK cutoff j2^|. However, we can rule it out if 
we make the additional assumption that the acceleration mechanism for cosmic rays is by 
"ordinary" electromagnetic fields. Since the glueball is neutral, it could never achieve more 
energy than the proton (or other charged hadron) which radiated it, and so the energy limit 
for the proton is the energy limit for any hadron, within this assumption and assuming as 
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well that no other charged hadron would remain stable. This seems safe, as all charged 
hadrons have valence quarks which give them quark content of order or larger than that of 
the proton. 

In the opposite case, e' < 0, where gluons propagate faster, a quark-rich hadron is ener- 
getically favored. The pion may be, and the kaon is expected to be, slightly quark-richer 
than the proton, but we find that neither is richer by enough to compete with the process 
p — > pD°, with the D° a cu or uc bound state. The kinematic limit for the process, from 
Eq. (|5.5|) . allows it to occur whenever e' < —1.8 x 10 -21 . The analogous process with a 
-D ± requires 6 MeV more energy in the center of mass frame, because of m D + — m D o and 
m n — m p . We have crudely estimated that, when 6 MeV is available in the center of mass 
frame, the center-of-mass lifetime for the process p — > pD° is around 10~ 9 seconds. The 
p — > pD° process occurs before there is enough energy to permit p —>■ nD + , provided the 
cosmic ray acceleration takes less than a few hundred seconds (in our frame), which seems 
reasonable for most proposed cosmic ray sources. Therefore the proton will lose energy to a 
neutral particle, which will not be accelerated further. The maximum energy a hadron can 
obtain is therefore the energy of a proton at the kinematic limit for this reaction. This leads 
to the limits on e' given in the introduction. 

We already saw that we can derive a limit when c 7 < (c q + c g )/2. We can also set a limit 
in the opposite follows. The proton has a larger photon (electroweak boson) content 

than the neutron, and when c 7 > c q ,c g , this means the limiting speed of a proton is faster 
than that of a neutron, opening up the process p —>■ ne + v e at high enough energies. (The 
reason the proton has more gauge boson than the neutron is, that it has more up quarks, 
and the electric charge of the up quark is larger than that of the down quark; so the quarks 
in a proton radiate more electroweak bosons.) The difference in electroweak boson content, 
barely discernible in Table HI is 0.14% of the proton's momentum. Neglecting for a moment 
the electron and neutrino masses, this is enough to counter the neutron's larger mass when 
e > .0014m p (m n — m p )/p 2 . Naively, this leads to a bound of 2 x 10 -23 . 

However, this is not the kinematic limit we find, because the timescale for the process 
p — > ne + v e involves the fifth power of the available energy. Together with the time dilation 
of the reaction, this leads to a time scale for the reaction, when close to the kinematic limit, 
of order 10 14 seconds. This is shorter than the likely propagation time of a cosmic ray, 
but it is much longer than the expected acceleration time of a cosmic ray, which we will 
take to be 10 4 seconds (of order the time variability of an active galactic nucleus). This 
knocks two orders of magnitude off our limit. If the correct time scale is shorter, our limit 
deteriorates as the 1/5 power of the inverse acceleration time. Note also that this limit relies 
on our assumption of ordinary acceleration. However, the energy involved is enough that 
the electron mass is negligible, and corrections from different leptonic dispersion relations 
will be too small to matter, according to constraints we will derive next. 

We have already seen how p — ► pe + e~ gives a strong constraint on the speed difference 
of g + q and leptons, when the leptons are slower. Similar considerations for p — > pvv prove 
that this limit applies both to singlet and doublet neutrinos. However, the limits on the case 
where the leptons propagate more slowly are much weaker; the analog of the constraint we 
found for the photon is 15 times weaker, because the lepton content of a hadron is second 
order in a, so the difference between the leptonic content of a neutron and proton is 10~ 4 . 
We actually get a stronger bound by combining the two-sided constraints on c 7 — (c q + c g )/2, 
just derived, with the Crab constraints 0] on the electron-photon velocity difference. This 



15 



gives the constraints on the c e , c p differences quoted in the introduction. 

The case where different quarks have different limiting velocities is richer, and may be 
phenomenologically viable; for instance, if the up quark is slower than all other particles, 
p — > A ++ 7r _ may be allowed, with most of the energy carried by the A ++ . If, on the other 
hand, the slowest propagating quark is one of the heavier species, say, the charm quark, 
then p — > pD° again becomes favored and the proton loses energy to neutral particles; so 
comparable bounds would exist in this case. 



B. Dimension 5 operators 

Next, consider the case with nonvanishing dimension 5 operators. Since they are odd 
under CPT, these operators always imply opposite dispersion corrections for a particle and 
its CP conjugate, so if the right-handed (SU(2) singlet) up quark has energy E UR = p+k UR p 2 , 
then the left-handed (SU(2) singlet) anti-up has energy E Ur = p — k UR p 2 . Similarly, opposite 
polarizations of gauge bosons have opposite corrections. The coefficient k is independent for 
each distinct Standard Model field; in particular, it is different for right and left handed up 
quarks, though it is the same for left-handed up and down quarks. 

Since we consider this possibility rather unlikely, we will not attempt to be very general 
in our considerations. Rather, we will assume that the order of magnitude of k is the same 
for each species, or at least the same within species of a given spin, and try to set some 
limits on it. 

We do not know the spin dependent PDFs of the proton very accurately. What we do 
know is that the x 2 moments of PDFs in hadrons are generally quite small, < 0.05, because 
most of a hadron's structure at the scale 10 11 GeV lies in very small x partons. Therefore, 
if k is comparable for different Standard Model species, then only the dispersion corrections 
of leptons and electroweak bosons will be important. 

With this in mind, consider again the process p —> pj. Take the polarization dependent 
energy shift for the photon to be = p ± k^p 2 , with ± according to whether the photon 
is left or right circularly polarized. The lower energy photon can always be produced; if 
the helicity of the proton forbids its production, we may consider p — > pyy, with an extra, 
low energy photon of the opposite polarization. The energy difference for collinear emission, 
with a photon momentum fraction of x, is 

Tn 2 Tn 2 
AE = p + — - — (1 — x)p - - t- p — r - xp+ \kJx 2 p 2 , 
2p 2p(l — xj 

2 

= _HhL^ + \ k \ x y . (5. 6 ) 

2p 1 — x 

This first vanishes when \k\ = 2m 2 /p 3 , x = 1/2. Writing k = 1/A, the scale A is p 3 /2m 2 ~ 
1.5 x 10 34 . The result should be slightly less than this, because the structure of the photon 
is not purely an x = 1 photon; but looking at Table ITT1 we see that not much of the structure 
has cascaded by this scale, so the correction is quite modest. 

We arrive at essentially the same constraint on left-handed leptons by considering p — > 
pvv\ this applies in each generation, since all neutrinos are light. For right handed leptons, 
the constraints involving electrons are the same, and for muons are weaker only by a factor 
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of 2/3, but for r leptons they are weaker by a factor of about 20, due to the mass of the 
particle to be produced. 

In the case where dispersion corrections are more than 30 times larger for QCD degrees 
of freedom than for leptons or electroweak bosons, these constraints do not apply. New 
constraints may or may not emerge, depending on exactly what the corrections are. For 
instance, if k g i uori is much the largest, then p — > Agg, with gg a spin-2 glueball, may be 
favored. If on the other hand the right handed up quark has the largest k, then p — > A ++ 7r - 
will probably occur, and this possibility cannot be excluded on the basis of cosmic rays, 
since the A ++ is charged and would be the most stable particle at large momenta. 

C. Dimension 6 operators 

If Lorentz symmetry is violated at a high scale, it is difficult to see how to avoid a 
dispersion correction E ~ p + 0(p 3 /A 2 ), with A a fundamental scale. Presumably A < m p \ 
the Planck scale, m p i ~ 1.2 x 10 19 GeV. For instance, if spacetime is a lattice, then lattice 
dispersion relations of this form are generically generated for free fields. 4 

Simplemindedly, one might expect the most natural form for dimension 6 dispersion 
corrections to be a universal K\ for all gauge fields and a Ki/ 2 for all fermions. Let us see 
what limits can be set in this case. If they are of comparable size, then a negative K\ leads 
to photon production, p — > pj, while a negative K 1 / 2 leads to p — > pe + e~ (or neutrinos). 
Table |H] shows that only 10% of a photon's p 3 weighted content is fermionic, while for the 
electron less than 5% is gauge boson; while both particles have over 100 times the size of 
dimension 6 correction that hadrons have. We therefore neglect the "cross-contamination" 
of 7 and e, and drop the small dimension 6 content of the hadrons. For the process p — > py, 
the process occurs when the energy shift, 

2 

AE = UhL^L + o.7Kx 3 p 3 , (5.7) 
2p 1 — x 

is positive. Here 0.7 is J x 3, y(x)dx for a photon, from Table |H] The resulting bound is, 
Ki > —1/ (4 x 10 22 GeV) -2 . For the case p — > pe + e~ , we may neglect the electron mass and 
consider the case where one electron emerges with very small x; approximately the same 
limit is obtained, K\/2 > — 1/(4 x 10 22 GeV) -2 . The scale appearing here is more than 2 
orders of magnitude larger than the Planck scale. 

For the opposite case, K\ or K 1 / 2 larger than 0, the proton energy rises faster than linear 
with p at large momenta, so ordinary Cherenkov radiation of soft photons is allowed, p — > py, 
above a threshold of 

2 

Tfl 

.0008^ + .0047ir 1/2 = — , (5.8) 

op 4 

where the numerical constants are the gauge boson and fermion dimension 6 content of the 
proton from Table ITT1 The numerical values of the bounds are, K\ < (6 x 10 21 GeV) -2 , 



4 Actually, if spacetime is a lattice, then interactions will generically induce O(a) corrections in dimension 4 
Lorentz violating operators. We have already seen that such corrections are ruled out by about 20 orders 
of magnitude. 
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K1/2 < (1.6 x 10 22 GeV) 2 . Surprisingly, these bounds are only about 4 times weaker; this 
is due to the more favorable kinematics for the p — > pj reaction in the positive K scenario. 

What if the dimension 6 operators show more species dependence? The cases < 0, 
-ft'iepton < have the same constraints just found; so do the cases K g > 0, K q > 0. For the 
case K q < but still universal between flavors, the dominant reaction is p — > pD°, which 
is allowed by virtue of the larger x quark content of the D° meson. The bound, derived 
in perfect analogy with the bounds found so far, is K q > — (6 x 10 21 GeV) -2 . If K only 
applies to first generation quarks, the strongest bound is from p — ► p7r°, and the bound 
is about 3 times weaker (but still exists). A similar bound on K g arises from p — > pD°, 
Kg > —(2 x 10 21 GeV) -2 . If we assume K q , K g = 0, we can also get bounds on positive K 1} 
K e from p — > p-f (with 7 soft), because of the photon and electron content of the proton; 
these bounds are K n < (3 x 10 21 GeV)" 2 , K e < (1 x 10 21 GeV)" 2 . If K g , K q are nonzero, 
then the bounds are, K 7 < 25\K q \, K e < 500|Kq| (or the above bound, whichever is weaker). 
Even these much weaker bounds suggest a trans-Planckian scale. 

VI. CONCLUSIONS 

We have shown that the mere existence of hadronic, high energy cosmic rays places 
constraints on Lorentz violating extensions of the Standard Model which are tighter than 
any previous constraints for a number of operators. To find these constraints, we first argued 
that Lorentz violating operators modify particle dispersion relations in a way determined by 
the partonic content of the particle in question. Then, we determined the partonic content of 
various hadrons, the photon, and the electron at the relevant scale, set by the energy of the 
cosmic rays in question. We then looked for Cherenkov-like processes made kinematically 
possible by these dispersion corrections, and excluded Lorentz violating processes which 
would prevent extremely high energy, hadronic cosmic rays from existing. 

We find that, generally, species dependence in the maximal propagation speed of a par- 
ticle, induced by dimension 4 Lorentz violating operators, is constrained to be smaller than 
a part in 10 21 . This applied for surprisingly many such dispersion corrections, and all of our 
limits are two-sided. Limits on dimension 5, CPT violating operators imply that, if gener- 
ated with 0(1) coefficients at the scale of Lorentz symmetry breaking, their intrinsic scale 
must exceed 10 34 GeV; for dimension 6 operators, the intrinsic scale must exceed 10 21 GeV. 
It is easy to evade the bounds on dimension 5 operators; they are absent in theories with 
CPT invariance. The limits on propagation speed differences and dimension 6 operators are 
harder to avoid in a theory with maximal Lorentz violation in the UV. (Similarly, maximum 
speed differences in theories without Lorentz symmetry are expected to be 0(a) ~ 10 -2 , 
as recently emphasized by Collins et. al. [IH, but it might be possible that there is some 
loophole in this argument. The argument does not apply if Lorentz symmetry is violated 
in some hidden sector, in which case the size of Lorentz violating operators in the Standard 
Model could be much smaller.) 

We should finish by commenting on some loopholes and limitations of our analysis. First, 
the loopholes. Most but not all combinations of Lorentz violating operators can be excluded. 
In particular, we find the interesting possibility that, if right-handed up quarks have a slower 
propagation speed than any other particle (either due to dimension 4 corrections or higher 
order in momentum dependent corrections), then the A ++ particle can be absolutely stable 
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at high energies. This possibility may also evade the GZK cutoff. This is an example of 
the narrow, special corners of the parameter space of Lorentz violating effects which are 
consistent with observed high energy cosmic ray events. 

Also note that not all forms of Lorentz violation lead to dispersion corrections of the 
form that we have bounded; for instance, spacetime noncommutativity can be understood 
in terms of the introduction of high dimension, Lorentz violating operators j^], but at 
leading order, the operators in question have no twist-two content (see [13)) so they do not 
modify dispersion relations. We therefore do not obtain strong constraints on spacetime 
noncommutativity from our results (although this will probably change if we work beyond 
leading order in the parameter describing the noncommutativity). 

One major limitation of our work is that we have considered only the case where there is 
a single preferred direction; Kostelecky's Standard Model Extension [l| allows more general 
corrections. Second, we have not been very exhaustive about what sorts of Lorentz violating 
operators to allow, for instance, generally taking the corrections of all quark flavors to be 
the same. Third, we have worked within the Standard Model, neglecting the possibility, for 
instance, that supersymmetry is realized at an energy scale not far above the electroweak one. 
This latter possibility would require the inclusion of several additional parton distribution 
functions in our evolution; we would find that the proton and other hadrons have nontrivial 
gluino and squark content. The tendency of hadronic PDFs to resemble each other would also 
be strengthened, because a s would run more slowly and the set of splitting processes available 
would be larger; therefore limits on, for instance, the quark-gluon velocity difference would 
be weakened. However, in general we do not believe that this case would be very different 
from the case of the Standard Model, because we expect the Lorentz violating operators to 
respect supersymmetry-otherwise, how would we obtain low energy supersymmetry in the 
first place? We leave these issues to future work. 
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APPENDIX A: SPLITTING FUNCTIONS AND COEFFICIENTS IN THE 
STANDARD MODEL 

The general form of the DGLAP equations is, 



Here fi(x,q 2 ) are the parton distribution functions; j is an index over distinct species. 
Pi<-j(y) is the ra t e at which species type j splits to produce species i with momentum 
fraction y. It is generally determined in a series expansion in the coupling, beginning at 




(Al) 



0(a/4n). 
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The DGLAP splitting functions, at leading order, are determined by the squared matrix 
element for a particle of momentum p to go to a particle of momentum xp + pj_ and a 
particle of momentum (1 — x)p — p±, evaluated at leading order in p\. Namely, if the matrix 
element (using the on-shell spinors or perpendicular polarization states as if the particles 
were on shell) squared, summed on final and averaged over initial spin states, is A4 2 , then the 
splitting function is given by M 2 p 3 x(l — x)/2[p x (rrp+p^)] 2 . That is, the matrix element 
is to be divided by the squared invariant describing how non-collinear the final states are, 
and multiplied by the energies of the incoming and outgoing states (when using relativistic 
state normalization). 

For a Lagrangian containing a gauge field and N{ Dirac fermions in representation / 
under this field, 

1 

we evaluate the splitting function for the process ip — > ^7, with 7 carrying momentum 
fraction x, by evaluating 

2 

u[(l-x)p + p±,af]tu\p,<Ti] , (A3) 

averaging over initial spins and colors and summing over final spins and colors. Here u\p, <jj] 
is the spinor for a particle of momentum p and spin a. A straightforward and well known 
calculation leads to the splitting functions, 



L = —F, U F^ + i>(fl + rr# , (A2) 



P^(y) = 2C F (j^^ - 1 - v + - y) 

P^(y) = U F N f (l-2y + 2y 2 ), 
P^(y) = 2C F Q-2 + ^ , 



P^(y) = AC A ) + \ ~ 2 + V - /) + (y C A - J F N^J 5(1 - y) , (A4) 

each times a/ 'Air. Here C F is the quadratic Casimir of the fermionic representation (C F = 4/3 
for SU(3), 3/4 for SU(2), and q 2 for U(l)), Ca is the Casimir of the adjoint representation 
(A^ c for SU(A C )), and t F is the fermion trace normalization (1/2 for SU(A C ), q 2 for U(l)). 
Contrary to common usage, the PDF for the fermion is taken in the above to be the sum of 
the usual fermion and anti-fermion PDF. 

If instead we consider a complex scalar field, the quantity to compute is 

- 2 

(p + xp + p±) ■ e , (A5) 
and splitting functions turn out to be a/4n times 
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The DGLAP splitting functions arising from the Lagrangian terms, 



L = {D^H)\D»H) + y t {Ht L P R t R + h.c.) , (A7) 

arise from evaluating \u(xp + Pi_)u{p)\ 2 . For the case where all three fields, H, t R , and 
are singlets, the splitting functions are y 2 /16n 2 times 



1 



P^(y) = (1 -</)(!-<%)- -5(1 



P^i(y) = y, 

P^(y) = -5(1 - y) , (A8) 

where the indices i, j refer to whether the field is t R or t^. 

The structure of Eq. (jAl|) means that the integral moments of interest to us obey closed 
sets of ordinary differential RG equations. This point is "well known to the people who 
know it well," but may not be familiar to all readers, so we go through it in a little detail. 
Introduce ^ 

fUq 2 )^ [ x n f t {x,q 2 )dx, (A9) 
Jo 

the moment of the PDF corresponding to the expectation value of a dimension ri + 3, twist-2 
operator. We obtain an evolution equation for /j jn from Eq. (jAlj) by integrating it over x n dx; 

/ x — dx = > / x dx —Pi^jiy) fj{x/y,q ) 

Jo a( i j Jo Jx y 

« 2 4M9 2 ) = W'* F X "dx Pl ^(y) fi(*/y,1 2 ) 

■ Jo y Jo 



dq 2 



fy n dy I ' z n dzP^ 3 (y)f 3 {z,q 2 ) 
j Jo Jo 

^j\ n j)(z,q 2 )dz^j ^ y n dyP^(y)dy^ 

y~\Cij,nfj,n{q 2 ) , with C ijin = / y n P i+ 

Jo 



i{y)dy. (AlO) 

The evolution equations for the fi^ n form a closed set of ordinary differential renormalization 
group equations, which mix particle species but not operator dimension n. 

It should not be surprising that the expectation values of all dimension D, twist 2 op- 
erators should satisfy closed renormalization group equations in this way. We expect an 
operator to mix under renormalization group flow only with operators of the same spin and 
equal or lower dimension. Since twist 2 operators have the maximum spin available at a 
given dimension, the twist 2 operators of dimension D can only mix with other twist 2 
operators at dimension D. 

It remains to evaluate the CV,- for all particles of the Standard Model, from the splitting 
functions we have just given. This is a straightforward combinatorial problem. We will write 
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the fermionic fields as [E, L, U, D, Q] for the right-handed electron, leptonic doublet, right- 
handed up, right-handed down, and quark doublet fields, [B, W, G] for the hypercharge, 
weak SU(2), and strong SU(3) (gluon) gauge bosons, and H for the Higgs doublet. The 
couplings are [a', a, a s , a y ], with a = a/4n and a y = y 2 /16n 2 with y the Yukawa coupling of 
the top quark, mf = y 2 v 2 /2 (v the Higgs vacuum expectation value, G F = l/v 2 \/2). The 
evolution equations for the integral moments of the PDFs are, at leading order, 



<fd 
dq 

q 2 d T (I , 3 \ _ 1 __. 1 



-y^E m — a'CffE m + a'CfgB , 



^ Lm = ^_ a > + _ a j CffLm + -aC fg W + -a'C fg B , 
dq 2 

dq 2 m 



D m = (^a s + ^C ff D m + ^C fg a s G + ^C fg a'B, 
{t as + t a ') ° ffUm + \ C f° asG + t° fg a ' B 



q 2 d 



dq 

fd 
dq 2 



+a y S m3 [2Cff y U 3 + Cf 2 f y Q3 + Cf s H] , 
2 Q m = (^a s + ^a + ^a'^j C ff Q m + ^C fg a s G + ^C fg aW + ^C fg a'B 



+0>ySm3 [CffyQ?, + 2Cf2f y U 3 + Cf s H] , 

B = a'C g f ^2 yEm + -^L rn + -D m + -U m + T^Qm j H ^ a 'C gg fB 



+-a'C gs H + -a'C ggs B , 



q 2 d Tjr ^ ^ f3 r ,3 
Jq 2 

+-aC gs H + -aC gga W. 



W = aC gf J2 {^L m + \q^\ + 2N g aC ggf W + 2aC ggg W 



dq 2 



^ (A 4 4 \ 

= a s C g f ( -^U m + - _D m + -Q m J + 2N g a s C gg fG + 3a s C ggg G , 



0/7 = (j a + \ a ') C ss 9 H + \aC sg W + l -a'C sg B 

+a y [C ssy H + C sf (Q 3 + 2U 3 )} , (All) 

where m is a generation index, iV s = 3 is the number of generations, and the C xxx coefficients 
are integrals of splitting functions given explicitly by, 

g 2 2 

g ' = - , Cggf — ) Cggg = , 

9 c — — r Zll 

111 

= g , C s / = - , Cf s = - , C ssy = -1 , (A12) 
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for the case of the dimension 4 integral moments and 
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Cf s = 


i C -- 


-1, 



for the case of dimension 6 integral moments. The a 2 s corrections can be found in [18[, and 
the Standard Model beta functions are in [26], and we will not repeat them here. 
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